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Abstract 

The epitaxial growth of {111} oriented Au on M 0 S 2 is well documented despite the large lattice 
mismatch (~ 8% biaxial strain), and the fact that a Au {001} orientation results in much less elastic 
strain. An analysis based on density functional and linear elasticity theories reveals that the {111} 
orientation is stabilized by a combination of favorable surface and interfacial contributions to the 
energy, and the compliance of the first layer of the MoS 2 . 
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FIG. 1: (color online) (a) The strain for An (left) and M 0 S 2 (right) in {111} orientation. The 
Molybdenum atoms are represented by purple circles, the Sulfur atoms by yellow circles and the 
Gold atoms by orange circles. a(Au) and a(MoS 2 ) are the lattice constants in Au (111) plane and 
M 0 S 2 respectively. The arrows indicate the biaxial strain in this case, (b) The strain in {100} 
orientation. The subscripts ”1” and ”2” are two perpendicular directions in the plane. ai and a 2 
are the lattice constants in corresponding directions. In this case, the sign and amount of strain 
in the two directions is different. 


The electronic and optical properties of transition metal dichalcogenides show much 
promise for technological applications [l|. Incorporating this material within devices will 
require either the growth of the dichalcogenides on other substrates, or growth of other 
materials on a dichalcogenide substrate. In this respect, the growth of Au on M 0 S 2 can be 
viewed as a prototypical system. 


The growth of Au on M 0 S 2 was studiec 
situ transmission electron microscopy 


in the mid to late 1960’s using early in situ and ex 
. Au was deposited on M 0 S 2 using evaporation, 
and was discovered to grow predominantly with a plate-like geometry in a {111} orientation. 
The (110) directions of the Au platelets were nearly aligned aloM the (1120) direction of 
the substrates. Similar orientations were observed for Ag nuclei 6|. 

The gro wth and evolution of the Au dims were modeled using theory available at the 
time j^, 7]. This theory, however, was rooted in an phenomenological understanding of 
the interfacial properties. Moreover, the substrate was treated as a typical bulk, and no 
accounting for the influence of the van der 'VVaals (VD’W) bonding within the substrate 
layers was attempted. Further, the Au clusters were approximated as spherical caps, as the 
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FIG. 2: (color online) The strain energy per unit cell for Au deposition on M 0 S 2 in two orientations 
computed assuming that the M 0 S 2 is unrelaxed, and that surface and interfacial energies make 
negligible contributions to the energy. Note that the strain energy associated with the {111} 
orientation of Au exceeds the strain energy of the {001} oriented growth by a substantial amount. 

TEM images did not allow for measurement of island thicknesses. Finally, only the {111} 
orientation was considered in any detail. 

In the most simple model of epitaxy, the substrate is assumed inhnite, and as a con¬ 
sequence, it does not relax during the growth process. For Au on M 0 S 2 , two possible 
orientations of the growing him are shown in Fig. [H For the {111} orientation, the biaxial 
strain required in the him is approximately 8%. In contrast, the {001} oriented him the 
strains are approximately -6% and 8% in the directions shown. Based on these strains, one 
would expect that the {001} orientation would be much more favorable, and consequently 
would be predominantly the experimentally observed orientation. Figure [2] compares the 
elastic energies of the two hlms (neglecting surface and interfacial stresses and energies) as 
a function of the number of layers of Au grown assuming the substrate is rigid. Clearly, 
the elastic energy of the {111} orientation is much larger than that of the {001} orienta¬ 
tion. This observation raises a fundamental questions regarding the growth: Why is the 
predominant orientation {111} and not {001}? 

In the following paragraphs a model to explain the experimentally observed him orienta¬ 
tion is developed. The model assumes that due to the weak VDW bonding between layers 
of M 0 S 2 , the surface layer is able to relax nearly independently of the remaining bulk lay¬ 
ers. The compliance of the substrate, when coupled with the surface and interface energies 
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(including strain energies) results in a lower formation energy for the {111} orientation as 
compared with the {001} orientation, in agreement with experiment. 

The model is built on two types of calculations. First, a continuum linear elastic model 
for the epitaxial growth of An on M 0 S 2 including the relaxation of the substrate and the 
influence of surface/interface energies and strains is developed. The parameters for this 
model are then determined using density functional theory (DFT) based electronic structure 
total energy methods. 

The continuum linear elastic model is developed by proposing a “synthesis path” and then 
computing the energy contributions along this path. (Note that this path is not necessarily 
experimentally accessible. It simply facilitates computation.) The initial step in forming an 
epitaxial layer can be taken as the extraction of a thin slab of the growing material from a 
bulk crystal. This extraction creates two surfaces, and the energy of the slab is increased by 
the surface energies (that reflect any strain in the as produced surfaces). The next step is to 
strain the him to its hnal strain in its epitaxially bonded state. This contributes the strain 
energy of the bulk plus any contribution to the strain energy from the surfaces. The third 
step is to separate the hrst substrate layer from rest of the substrate to create a freestanding 
M 0 S 2 layer. This adds a layer separation energy to the system. In the fourth step, this 
free-standing layer is then strained introducing the strain energy of the single layer. In the 
hfth step, the epitaxial him is welded to the free standing layer. This has the consequence 
that one of the strained An surfaces is replaced with a strained AU-M 0 S 2 interface. In the 
hnal step, the hlm/free-standing layer is readhered to the substrate, returning the layer 
separation energy to the thermal bath, but introducing the energy required to slip the top 
M 0 S 2 layer relative to the remaining layers. Taken in total, the sum of the changes in energy 
for both the substrate and the epitaxial him as compared with their bulk counterparts, AE 
can be written: 

AE T EfiiYn ~\~ Egy^i) EIg T Eglip (1) 

where Eau^sut is the energy of the strained An surface, Efu^ is the energy of the strained 
him neglecting surface contributions, Egut is the strain energy of the hrst substrate layer, 
EAu/s,int interfacial and strain energy of the AU-M 0 S 2 substrate interface, and EgUp is 

the slip energy between the hrst layer of the substrate and the remaining substrate. 

The hrst four contributions to AE can be expressed analytically using linear elasticity 
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theory: 



1 

( 2 ) 

E film 


Esub 

_ 1 

2^^sub,ijkl^sub,ij^sub,kl-^sub 

(3) 

Eau^suv 


(4) 

EAuj S^in 

= 7r + + 2^I,ijkl^I,ij^I,klAl 

(5) 


The Cf^ijki, Csub,ijki are the elastic constants for the him and the top layer substrate. The 
strain tensors are indexed similarly. We approximate the surface/interface stress energies up 
to the second order of the strain tensor: 7 ^ and 7 / are the unstrained surface and interfacial 
energies, respectively; fs,ij and fi^ij are the linear surface stress terms while Cs,ijki and Cj^ijki 
are the quadratic terms or the effective elastic constants of the surface/interface. The top 
layer of the substrate is treated as a 2D material. (We assume that the VDW interaction 
within the substrate is sufficient to insure the substrate remains hat during the epitaxial 
growth). Dehne Vf to be the equilibrium volume the him would have if it were part of a 
bulk An crystal. Aguh, Ag and Aj are the reference unit cell areas for the substrate, interface 
and surface strains respectively. Ag and Aj are taken equal to Af, the area covered by the 
An him with its bulk lattice parameter. Agub is taken to be the equilibrium area of the 
monolayer M 0 S 2 unit cell. 

The slip energy arises from displacing the top substrate layer relative to layers below. 
An approximation for this energy is made by investigating a similar slip in a bilayer M 0 S 2 
system. Two M 0 S 2 layers are placed relative to each other in the same way as the two 
adjacent layers in the bulk M 0 S 2 . The slip system is then created by straining one layer by 
5% biaxially while keeping the other one hxed. (5% is about the typical amount of strain in 
the AU-M 0 S 2 epitaxy system.) The interlayer distance in both cases is hxed at 6.25 A , the 
value obtained from DFT calculations of the pristine bilayer M 0 S 2 system. The slip energy 
consists of only the change in the interlayer VDW energy since the change of the VDW 
energy within the strained layer is captured in the elastic energy computation. The VDW 
interaction energy is directly calculated following Grimme’s D2 method 8 |. The diherence 
of the interlayer VDW energies between the slipped system and original system gives an 
estimate for the slip energy. The computed slip energy in this case is less than 2 meV per 
unit cell and will make a negligible contribution to AE. (This assertion is borne out by 
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the calculations below. For more details of the slip energy calculation please refer to the 
Supplemental Material.) Therefore, for simplicity, this term is neglected in the model. 

The value of AE at equilibrium is then determined by minimizing the right hand side 
of Eqn. dl]) with respect to the strains in the him and the hrst substrate layer subject to 
the constraint that the hrst substrate layer and the him are lattice matched across their 
interface. (A more detailed description is provided in the Supplemental Material.) 

The parameters that enter the theory are determined either by htting the results of DFT 
based electronic structure total energy calculations to the continuum theory or by direct 
DFT calculations. DFT calculations are performed using the plane-wave code VASP {q]. 
The exchange and correlation energy are described by generalized gradient approximation 


proposed by Perdew, Burke, and Ernzerhof 


Q, 


Electron-ion interactions are treated with 


projector augmented wave potentials [n|. All calculations are performed using a plane-wave 
basis with a 350 eV energy cutoh. The precision tag is set to “accurate.” The convergence 
criterion for self-consistent held loop is 1 x 10“® eV. A 20 A vacuum slab is added along the 
direction normal to the growth plane to separate the system from its periodic image. 

Two epitaxial conhgurations are considered, as shown in Fig. |3l In {111} epitaxy 
(Fig. [3](a) and top right of [3](c)), An (111) plane is lattice matched onto the M 0 S 2 . In 
Fig-EKc), the green circles are the gold atoms in the hrst layer that he exactly on top of the 
Sulfur atoms. The green dashed lines dehne the unit cell. The gold atoms in the second 
and third layer are projected onto the growth plane and are represented by the blue and 
red circles, respectively. The stacking of gold layers follows the “ABCABC...” sequence to 
be consistent with the Face Center Cubic (FCC) structure, with the “B” sites residing on 
top of the Mo atoms. (The “ACBACB” sequence is slightly higher in energy (about 0.01 
eV/unit cell) than the “ABC” conhguration. Thus we take the “ABC” sequence for {111} 
epitaxy, though both structures have been realized experimentally jo].) In {001} epitaxy, 
the An (001) plane is lattice matched onto the M 0 S 2 (Fig. [3](b) and bottom left of |3](c)). 
In the top view, the gold atoms in the hrst layer are represented by grey circles. The grey 
dashed line dehnes the unit cell. The black circles are the projections of the gold atoms 
in the second layer. The stacking sequence follows “A'B'A'B'...” sequence. We constructed 
models for both conhgurations with the number of An layers vary from 3 to 30. A 14x14x1 
Monkhorst-Pack k-point grid is used to sample the Brillioun zone for {111} epitaxy. In 
{001} epitaxy, a 13x8x1 Monkhorst-Pack k-point grid is used. These structures are relaxed 
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FIG. 3: (color online) Side view of (a) {111} epitaxy unit cell with 3 layers of Au and (b) {001} 
epitaxy unit cell with 3 layers of Au. The vacuum layer is used to prevent the interaction of the 
system with its periodic image, (c) Top view of two types of the epitaxy for Au grown on MoS 2 . 
The unit cell of {111} epitaxy is defined by the green dashed line . The unit cell of {001} epitaxy 
is defined by the grey dashed line. 

until the maximum Hellmann-Feynman force on any atom is below 0.01 eV/A. 

The DFT calculated lattice structures and elastic constants of bulk Au and monolayer 
M 0 S 2 are in good agreement with previous studies (details in the Supplementary Materi¬ 
als). The calculated elastic constants are used for the parameters Cf^ijki and Csub,ijki- The 
calculated lattice structures are used to determine V/, Agub, Ag and Aj. 

Within the density functional theory approach, the change in energy associated with our 
synthesis path can be written (neglecting the slip energy): 

AF/ Eiq^ IT'AuEh^Au E}jA1oS2i (6) 

with Etot the total energy of the lattice matched slab including one M 0 S 2 layer and uau 
layers of Au (in each unit cell). E^^au is the energy per layer of bulk unstrained Au and 
Eb,MoS 2 is the energy of the monolayer M 0 S 2 . We £t the results of DFT calculations to both 
Eqn. ([1]) and the equilibrium strains (with more weight on the energy £t side) to determine 
the parameters remaining in that expression. (These turn out to be linear combinations of 
the surface and interfacial energy terms. See the Supplemental Material for more details 
about htting processes.) 

Figure m compares the AE’s for {111} and {001} oriented epitaxial growth. The results 
of DFT calculations are displayed, as is their £t to Eqn. ([T]). It is apparent that the elasticity 
theory based model does an excellent job of describing the DFT results, and can be used to 
extrapolate behavior to thicker hlms. In addition. Fig. 0] displays the predictions assuming 
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a rigid substrate (traditional epitaxy), and the interfacial energies and stresses computed 
within the model. 

Examination of the Fig. 0] reveals two important points. First and foremost, within the 
compliant substrate epitaxial model, for all considered numbers of An layers, {111} epitaxy 
is energetically favored over {001} oriented epitaxy. However, the energy difference is not 
too large, and one might expect to see both orientations, with the {111} orientation favored. 
Thus the compliant substrate epitaxy model is able to explain the experimental observations 
of {111} oriented epitaxy of An on M 0 S 2 . 

Second, the variation of AF with thickness is sublinear in the case of compliant substrate 
epitaxy. This sublinear behavior originates in the fact that the strain energy per layer is 
decreasing as the him thickness increases. As the him gets thicker, it becomes elastically more 
stih, and the hrst layer of the substrate is forced to deform to a greater extent. Eventually, 
when the him is inhnitely thick, only the hrst substrate layer deforms, and the elastic energy 
saturates at a constant. To our knowledge, this reduction in strain energy of the him with 
increasing him thickness is a unique feature of compliant substrate epitaxy. 

These observations suggest the following understanding of the preferred orientation of the 
hlms. First, the relaxation of the substrate is not expected to change the sign of the strain 
energy diherence between {111} and {001} oriented hlms. Based on strain energy alone, one 
would still expect {001} oriented hlms, even for a compliant substrate. This implies that 
the sum of the {111} An surface and AU/M 0 S 2 interfacial energies is less than the sum of 
the {001} An surface and AU/M 0 S 2 interfacial energies in the compliant substrate epitaxy, a 
fact consistent with our htted parameters (See Supplemental Material). So in this instance, 
the properties of the interfaces and the surfaces dictate the orientation of the growing him. 

This observation is not, in itself surprising. When the hlms are nucleated, they are very 
thin, and the surfaces and interfaces can make a signihcant contribution to AE. However, 
the persistence of the orientation preference with increasing him thickness is remarkable. 
Typically, one expects that the strain energy diherence would become the dominant term 
in AE, and the favored orientation would change to {001} for thicker hlms. This is where 
the compliance of the substrate becomes important. Since the substrate is compliant, the 
strain energy of the An him no longer scales with the volume of the An him. Instead, 
it monotonically decreases with increasing him thickness, enabling the interfacial and sur¬ 
face contributions to the energy to dictate the stability of the two orientations for all him 



FIG. 4: (color online) A comparison between the elastic energy of compliant substrate epitaxy 
and that of traditional epitaxy that assumes a rigid substrate. Note that surface and interfacial 
energies and strain energies are included in the lines labelled traditional. DFT results are shown 
as dots, and the lines passing through them are the results of the continuum theory fitted to the 
DFT data. Inclusion of surface/interface effects alone does not stabilize the {111} orientation. 
This stabilization requires a compliant substrate. 


thicknesses. 

The model presented here is the not hrst proposed that exploits a compliant substrate. 
As early as 1991, Lo suggested that the quality of some epitaxial dims could be improved by 
employing a compliant substrate 12|. Such a substrate would enable relaxation of the him 
at the expense of the substrate, but had the potential to increase the quality of the epitaxial 
him. Lo suggested that such substrates could be produced using standard lithographic 
methods. Later, Jesser et al. proposed that a compliant substrate might be developed by 
introducing a subsurface twist boundary [^. Here, it is noted that layered materials with 
VDW bonding between them form naturally compliant substrates for epitaxial growth jl^ . 
Moreover, the VDW bonding, while enabling lateral slip of the substrate, will resist buckling 
of the him, and thereby help to improve its quality. Though we have examined An on M 0 S 2 
in detail, the idea is quite general, and should apply to other systems as well. 

Note that this version of compliant substrate epitaxy is diherent from the so-called van der 


Waals epitaxy 


15| . In the compliant substrate epitaxy model, lattice registry is maintained 


across the interface, and a degree of covalent bonding between the substrate and the him is 
allowed. In contrast, in the van der Waals epitaxy case, one does not expect lattice registry 
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between the film and the substrate. 


We note that compliant substrate epitaxy may have interesting implications for strain 
engineering and processing of thin dims. Consider the strain field near a small An island in 
the early stages of the him growth. The M 0 S 2 under the island will be strain in compression. 
This will naturally be accommodated by a tensile strain surrounding the island. Therefore, 
it might be possible to engineer the positions of the An nuclei so as to induce a desired 
strain pattern into the hrst layer of the M 0 S 2 substrate. Moreover, it might be possible to 
use the epitaxial binding as a means of freeing and manipulating the top layer of the M 0 S 2 , 


just as was demonstrated for C nanotubes 


16[ |. Finally, the fact that under conditions of 


compliant substrate epitaxy, surface and interfacial energies can dictate the orientations of 
thick dims raises the possibility of using surface chemistry and/or surfactants to engineer 
the relative stability of diderent dim orientations within the same materials system. 

In conclusion, the epitaxial growth of An on M 0 S 2 is studied using a combination of 
elasticity and density functional theories. It is shown that the compliance of the substrate 
in conjunction with the relative surface and interfacial energies and strain energies stabilize 
the {111} growth orientation despite the large lattice mismatch. This system is thus an 
example of compliant substrate epitaxy. 

This work was supported by the Director, Ofdce of Science, Ofdce of Basic Energy Sci¬ 
ences, Division of Materials Sciences and Engineering, of the U.S. Department of Energy 
under Contract No. DE-AC02-05CH11231. 
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1 . The lattice structures and elastic constants of Au and M0S2 from DFT calcula¬ 


tions. 


The lattice structures and elastic constants of bulk Au and monolayer M 0 S 2 computed 
using density functional theory are presented in this section. The parameters for the DFT 
calculations are same as those reported in the main text except where noted. For the 
bulk Au primitive unit cell, a symmetrized 12x12x12 Monkhorst-Pack k-point grid is used 
to sample the Brillioun zone. For the primitive cell of monolayer M 0 S 2 , a symmetrized 
14x14x1 Monkhorst-Pack k-point grid is used to sample the Brillioun zone. 

The computed Au lattice constant 4.168 A, and is in reasonable agreement with experi¬ 
ment result of 4.062 A and in good agreement with previous DFT result of 4.154 A [l|. The 
in-plane lattice constant of M 0 S 2 from our DFT calculation is 3.190 A and the thickness of 
the monolayer (the separation distance between the top and bottom Sulfur layers) is 3.130 
A. These are also in good agreement with the experimental values of 3.122 Aand 3.116 A, 
respectively j2|. The computed elastic constants for these two materials are summarized 
in Table [Tl Our calculated results for gold elastic constants (from GGA PBE) are softer 

n 

than the experimentally measured values [3|. We also calculated the Au’s elastic constants 
using PAW LDA potential (parameterized by Perdew and Zunger) { 4 ]. The LDA results 
overestimate the elastic constants of Au, as shown in Table H We further double checked 
our epitaxy system calculations (with 3, 21 and 24 layers of Au) using LDA. In the LDA 
calculations, the {111} epitaxy systems are still more stable than the {001} epitaxy systems, 
as shown in Fig. [H which is consistent with the GGA results. We, therefore, present GGA 
results in the main text. 

In the Ref [^, a discussion of the elastic constants of M 0 S 2 (both measured and predicted) 
is given. The G 12 measured experimentally is negative, which contrasts with the positive 
value computed here and with the DFT-D2 calculations done by Peelaers and van de Walle. 
Peelaers and van de Walle argued that the G 12 measurement is indirect, and might be 


complicated by the 


ubricating properties of M 0 S 2 . In order to directly compare our results 


with those in Ref js], we follow their convention of calculating the elastic constants for 2 D 
material, in which the thickness of monolayer M0S2 is chosen to be the half length of the 
unit cell vector perpendicular to the basal plane in bulk M0S2 (contains 2 M0S2 layers). 
Our calculated elastic constants are in good agreement with Peelaers and van de Walle’s 
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FIG. 1. (color online) IS.E calculated by LDA for the epitaxy systems with 3, 21 and 24 layers of 
An. The {111} epitaxy systems are lower in energy than the {001} epitaxy systems. 


TABLE I. The elastic constants of bulk An and M 0 S 2 . 



Elastic Constants 

GGA [GPall 

LDA [GPa] 

Experiment [GPa]^ 


Cii 

154 

219 

192 

Au 

C12 

136 

183 

163 


C44 

31 

41 

39 


Cii 

211 

237 

238 (238) 

M0S2 

C12 

53 

56 

-54 (64)1 


The GGA (PBE) results are used in the calculations described in 


the main text. 


The previous reports of Au elastic constants can be found in Ref 3| . The reported M 0 S 2 elastic 


cons 

Ref 


Ints is from Ref The DFT calculated elastic constnats from this reference is in parenthesis. 

I argued that the negative G 12 from the experiment is not measured directly and might come from 


lubricating properties of M 0 S 2 . 


results. The calculated structures and elastic constants are used for the parameters of our 
continuum elastic model, i. e. the equilibrium volume of the film, the equilibrium area of the 
substrate/surface/interface and the elastic constants of the film and the top layer substrate. 
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2 . The slip energy estimation in M0S2 substrate. 


In this section, an esimate of the hrst layer slip energy for the M 0 S 2 snbstrate is presented. 
The system we use is a bilayer M 0 S 2 system. A top view of this system is in Fig. [2] with 
only the lattice points shown. The two layers coincide at the origin point in the top view 
(left corner). Each layer contains 21x21 unit cells. The Mo and S atoms are added to the 
lattice points in the same way as that in the two adjacent layers in the bulk M 0 S 2 . The 
interlayer distance is 6.25 A, as obtained from DFT with van der Waals corrections of the 
pristine bilayer M 0 S 2 system. To create a slip similar to the compliant substrate epitaxy 
system, the hrst layer M 0 S 2 is strained biaxially by 5% while the second layer is kept hxed. 
This is the typical amount strain in the AU-M 0 S 2 epitaxy system. This strain results in 
a structure wherein the lattice parameters of the slipped layer and the unslipped layer are 
commensurate (the lattice points coincide in xy plane at the corners of the strained layer). 
When the strain is applied, the Mo and S atoms are also displaced with the lattice points 
and no internal relaxation of the unit cell is allowed. 


Biaxial strain on 
the black lattice 




FIG. 2. (color online) The two layer M 0 S 2 system (top view). Only the lattice points are shown. 
The black points indicate the first layer who has 5 % of biaxial strain. The blue points indicate 
the second layer which is without the strain. The lattice points of the two layers also coincide at 
the corners of the strained layer. 


We follows Grimme’s D2 method jo] to include the van der Waals (VDW) interaction 
energy and calculate the interlayer VDW interaction using the formula: 


Ecorr 


Natl Nat2 

-vv 


c, 


^6 fdmp{rij) 
ij 


( 1 ) 
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Where the index i goes over the atoms in the hrst layer and j goes over the atoms in the 
second layer. Tij is the distance between atom i and atom j. C^ij is a coefficient which 
depends on the types of atom i and atom j; 



( 2 ) 


The parameter C^i and Cqj are tabulated for each element and are insensitive to the par¬ 


ticular chemical situation. The values of Cg* for the elements in the hrst hve row of the 
periodic table are provided in Ref |6[. fdmpifij) is the damping function whose expression 
is: 



( 3 ) 


sg is a global scaling parameter which depends on DFT functionals used. We choose the 
optimized value for PBE functional (S'g = 0.75) since the structure of M 0 S 2 is determined 
by DFT calculations with PBE functional. The coefficient Roij also depends on the types 
of atom i and atom j: 


( 4 ) 



The values of Roi for the elements in the hrst hve row of the periodic table are also provided 
in Ref 6|. 


To estimate the slip energy, we compute the interlayer VDW energy of the original and 


slipped systems described above. The energy diherence between these two values provides 


an estimate of the slip energy. Only the interlayer interaction is counted since the change 
of the VDW interaction energy within the strained layer is captured in the elastic energy of 
the substrate. The VDW energy diherence calculated from the above method is 0.0019 eV 
per primitive unit cell of the M 0 S 2 . As mentioned in the main text, this slip energy makes 
a negligible contribution to AE. A similar order of magnitude of slip energy is expected in 
{001} epitaxy and other compliant substrate epitaxy systems. 

3. Detail strain analysis and minimization of AE in {111} and {001} compliant 
snbstrate epitaxy. 

As shown in Fig. 1(a) of the main text, the lattice mismatch between An and M 0 S 2 is 
about 8% biaxially in the case of {111} epitaxy. In this case, the An him has a biaxial tensile 
strain (e^^^n) in the growth plane and a normal strain component {eAu, 33 ) in the direction 
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perpendicular to the plane. The shear strain is zero. The strain of the interface and surface, 
as dehned in the main text, are same as the strain of the An him in the growth plane. The 
strain for the M0S2 is a biaxial compressive strain (eMoS2,ii) ^Iso with no shear component. 
With this form of the strain tensor given, the expression of AE can be written down based 
on Eqn. ( 1 ) to ( 5 ) of the main text. Furthermore, the value of eMoS2,ii is related to eAu,11 
through the condition of lattice matching at the interface: 

+ ^Au, 11 ) = flMoS2(l + CMoS2,1i) (5) 

Once the number of An layers is given, the only variables in the AF’s expression are eAu,ii 
and cam.ss. The AE can then be minimized with respect to eAu,ii eAu,‘iz- Moreover, in 
this case, only the combined linear term (/s + /r), the combined quadratic terms + 

Csub,1122 + C/piii + C/,1122) and the combined free interface and surface energies (7^ + 77-) 
enter into the minimized expression for AE. (Due to the symmetry of the surface/interface, 
the linear term can be taken as scalar j^.) Thus the combined terms can be htted as a 
whole to the DFT results. 

In the case of { 001 } epitaxy, the lattice mismatch is about 8% in direction 1 and -6% in 
direction 2, as shown in Fig. 1 (b) of the main text. In the growth plane, the An him has a 
pure tensile strain 6(4^ in direction 1 and a pure compression strain e(4„ 22 direction 2 . 
The An him also has a normal strain component e'Au 33 perpendicular to the growth plane. 
Similarly, the strain of the interface and surface are same as the strain in the An him growth 
plane. In the M0S2, likewise, it has normal compression strain 11 direction 1 and 
normal tensile strain 22 direction 2 . The shear components for both An and M0S2 
remain zero. There are now two constraint equations relating the strain of An and M0S2 in 
the growth plane: 

ai(An)(l + ^'au,Ii) ^ ®1(MoS 2)(1 + ^MoS2,ll ) (6) 

a2(An)(l + ^'au,22) = ®2(MoS'2)(l + ^'moS2,22 ) (7) 

Once the number of An layers is given, the variables in the total energy expressions are 
eV,22 and 6^,33- can then be minimized with respect to 6^,11, eV,22 and e^,33- 
In this case, we have one combined linear term (/^ + //), two combined elastic constants 
(C'sub,1111+07^1111 and 05.ufe^ii22+C'7_ii22) and the combined free surface and interface energies 
{is + 7/) entering the minimized AE expression. Just as in the case of { 111 } epitaxy, these 
parameters are htted as a whole to the DFT results. 
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In more detail, the fitting process includes both the energy and the in-plane strains 
through, These are incorporated within a sum of squares, 

y = w ^(ABli] - + (1 - y 5 ;^(£|j] - e\j])\ (8) 

* 3 

In the above equation, /S.E[i] is the DFT calculated excess energy of AU-M 0 S 2 system with 
i layers of An (dehned by Eqn. (6) of the main text). AA'[z] is the excess energy from the 
elasticity model also with i layers of An (the minimized Eqn. (1) of the main text). The 
strain sum term is dehned similarly. For {111} epitaxy, the in-plane biaxial strain comes 
into the strain sum. For {001} epitaxy, both en and €22 are included the strain sum. The 
expressions for the e[j] are also obtained from DFT results. The indices i and j run over all 
the computed DFT data, w is a weighting parameter that can be used to tune the relative 
importance of the ht to the energy versus the ht to the strain. This parameter is arbitrary. 
Since the DFT calculated energy is more accurate the DFT calculated strain, we choose to 
put more weight {w = 2/3) on the energy. The htting process minimizes with respect 
to these unknown parameters in the elasticity theory model. The htted results for both 
orientations are summarized in Table m 


TABLE II. The values of fitting parameters □ 



Linear terms 

Quadratic terms 

Free energy terms 


fs + fl 

CsubAlll + Csub,1122 + Cl nil + ( 7/1122 

li + ls 

{111} epitaxy 


0.167 

3.988 

0.067 


fs + fl 

CsubAlll + CsuhAl22 + (7/1122 

li + ls 

{001} epitaxy 


-0.541 

45.703 45.656 

0.099 


^ All fitting parameters are in unit of [eV/A^] 


The htted energy curves are shown in the Fig. 4 of the main text. The htted strain 
curves are shown in Fig. |3l The strain curves are also in reasonable agreement with the 
DFT predictions. 
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FIG. 3. (color online) The fitted An strain curve and DFT An strain data for (a) the biaxial strain 
in {111} epitaxy and (b) the two in-plane strains in {001} epitaxy. 
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